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Abstract 

We present a detailed description of the Q 2 and x dependence of 
the structure functions F2 and Fl as extracted from the Deep Inelastic 
Scattering data at HERA in the small Bjorken x region. Making use 
of a collinearly-improved BFKL equation at next-to-leading order and 
a treatment of the running of the coupling using non-Abelian physical 
renormalization together with the BLM scale choice allows us to reach 
low values of Q 2 . We also provide some predictions for future lepton- 
hadron colliders. 



1 Introduction & review of our framework 

Recently, in [lj, we calculated the effective slope A for the structure function 
F2 when parameterized as ~ x~ x at low Bjorken x, which ranges from A ~ 0.1 
to ~ 0.3 when moving from low to high values of Q 2 . The importance of a 
correct understanding of the HERA data using perturbative QCD techniques, 
in particular for the physics program at the Large Hadron Collider (LHC), has 
been largely discussed in the literature (see, e.g., Ref. [2]). 

In this Letter we will not only further investigate the Q 2 and x dependence 
of F 2 but also study in detail the longitudinal structure function Fl- As in 1 
we will explore the small x region using the next-to-leading order (NLO) 3 
BFKL [i] equation with collinear improvements, together with optimal renor 
malization schemes. 
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Let us first briefly review our formulae [Tj. At small x ~ Q 2 /s, with s being 
the squared center-of-mass energy, we can apply high energy factorization and 
write the structure functions Fj, I = 2, L as (note that the integrations take 
place in the transverse space with q = yqj) 

Fi(x,Q 2 ) = f^^( q ,Q 2 )^p( P ,Q 2 )^(x, q ,p). (1) 

J tt q 2 J T[p z v ' ' 

$p is the non-perturbative proton impact factor which we model using 

C I V 



where we have introduced two free parameters and a normalization. is the 
impact factor associated to the photon which in |l] we treated at leading-order 
(LO), i.e. 



nq 2 



where 



cAv) = — -. — — - — \tt sech(vrz^) tanh (in/) (4) 
4 (u + 



v = i(l/2 - 7), fi 2 = (11 + 12z/ 2 )/8, VL L = v 2 + 1/4, and the strong coupling 
a s is fixed at the renormalization scale fi 2 . In the present work we will also 
use the kinematically improved impact factors proposed in [5j|6], which include 
part of the higher order corrections by considering exact gluon kinematics. Its 
implementation requires to replace the functions ciiv) by 6/(7,0;) where 

~ ( 4r( 7 +e+i)r(i+ 7 ) [M7+O - tM (3^ 2 - e + 1) - 6^] 
Ci(7 ' w) - er(i + ^-5e + 4) (5) 



and c 2 = cl + Or, with 

r(7+0r(7) 



j-2£u; (£ 2 + 3 2 + 6w + 11) 



Cr(7 ' w) er(i+c)(e-5e+4) 

+ [V>(7 + - ^(7)] [e 4 - 10e 2 + 3w 2 (w 2 + 2u; + 4) - 2co (e 2 -l) + 9]J. (6) 

^(7) is the logarithmic derivative of the Euler Gamma function and £ = 1 — 
27 + u, while to is the Mellin variable conjugate to x in the definition of the 
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gluon Green function J 7 , see Eq. ^ below. The main difference between 
these impact factors is that the LO ones roughly double the value of their 
kinematically improved counterparts in the region with small \v\, while being 
very similar for \v\ > 1. 

The gluon Green function can be written in the form 

r( x , q , P ) = l f^L f^L ±( q -\\-» L_, (7) 

7r J 2m J 2iti q 2 \p 2 J oo - a s K, (7) ' 

with a s = a s N c /ir. The collinearly improved BFKL kernel as introduced in 
eq. is an operator consisting of a diagonal (scale invariant) piece x{l) with 
eigenvalue 

x{l) = ®sXo(l) + o?,X\(l) ~ 7,® 2 s Xo'(l)Xo(l) + Xrg(« s , 7, a, 6) , (8) 

where X o(l) = 2^(1) - Hi) ~ ^ ~ l), a = - - | and b = 

~ §iv ~ ~ T§> P ms a term Xrc( 7 ) proportional to j3o which contains the 
running coupling corrections of the NLO kernel [f\: 



Xrc(7) = ai^- [Xo(r/)d^- 9 7 xo(7) + 21og(^)j. (9) 

The precise form of the NLO kernel xi can be found in [TJ. The resummation 
of collinear logarithms of order and beyond is realized by the term [lj[8j[9] 

XRG(«s,7> a , b) = a s (l + aa s ) (^(7) - ^(7 - ba s )) 

C? 7T 2 1 °° / 

2 sm (tt 7 ) 2 ^ y 

2a s (l + aa. 



+ ^(7 - 1 - m + 6a s ) 2 + 4a s (l + aa s ) . (10) 
1 — 7 + m J 

Our final expression for the structure functions reads 

F,(x,QV/*^ [ 1 + g^|±^ l 0g (l) (11) 



x ( i(7rcoth(7rz/) — 27rtanh (7tz/) — M/(z/)) — ^ ( 8 — — — iv 



Ql 



2+ lv 

cAv) 



where M2 and Ml can be found in PH. For the kinematical improved version 
of Ff we replace Cj(v) by 6/(1/2 + zv, x(l/2 + ^)). 
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In Eq. (12) the scale of the running of the coupling has been set to 
= QQq. Building on the work of [lOj we found in [l] that in order to 



obtain a good description of the Q 2 dependence of the effective intercept of 
F 2 , A, for x < 10 -2 , it is very useful to operate with non-Abelian physical 
renormalization schemes using the Brodsky-Lepage-Mackenzie (BLM) optimal 
scale setting 11 with the momentum space (MOM) physical renormalization 
scheme 
reader to 




For technical details on our precise implementation we refer the 
see also 13 for a review on the subject and 14 for a related 



work). More qualitatively, in these schemes the pieces of the NLO BFKL ker- 
nel proportional to /3o are absorbed in a new definition of the running coupling 
in order to remove the infrared renormalon ambiguity Once this is done, the 
residual scheme dependence in this framework is very small. We also found 
it convenient |l] to introduce, in order to describe the data with small Q 2 , 
an analytic parametrization of the running coupling in the infrared proposed 
in (15|. 



2 Comparison to DIS experimental data 

In the following we compare our results with the experimental data for F 2 and 
F L . 



2.1 F 2 

Let us first compare the result obtained in [l] for the logarithmic derivative 



d\og F 2 /d \og(l/x) using Eq. (12) with a LO photon impact factor and our 
new calculation using the kinematically improved one. In Fig. [T] we present 
our results with the values of our best fits for both types of impact factors and 



compare them with the Hl-ZEUS combined data 16 for x < 10~ 2 . The values 
of the parameters defining the proton impact factor in ^ and the position 
of the (regularized) Landau pole (we use n/ = 4) for the strong coupling are 
5 = 8.4, Q = 0.28 GeV, A = 0.21 GeV for the LO order case and S = 6.5, 
Qo = 0.28 GeV, A = 0.21 GeV for the kinematically improved (note that the 
normalization C does not contribute to this quantity). 

The LO impact factor generates lower values than the kinematically im- 
proved one in the high Q 2 region and slightly higher ones when Q 2 < 2 GeV 2 . 
It is interesting to see how the approach presented here allows for a good de- 
scription of the data in a very wide range of Q 2 , not only for high values, where 
the experimental uncertainties are larger, but also in the non-perturbative re- 
gions due to our treatment of the running of the coupling. 
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Figure 1: Fit to A for F 2 with the LO photon impact factor (solid line) and the 
kinematically improved one (dashed line). The data set has been extracted 
from 16 



Encouraged by these positive results we now turn to investigate more dif- 
ferential distributions. We select data with fixed values of x and compare the 
Q 2 dependence of our theoretical predictions with them, now fixing the nor- 
malization for the LO impact factor to C — 1.50 and 2.39 for the kinematically 
improved. Our results are presented in Fig. [2} The equivalent comparison to 
data, this time fixing Q 2 and looking into the evolution in the x variable, is 
shown in Fig. [3| We observe that our predictions give a very accurate descrip- 
tion of the data for both types of impact factors. 

Let us remark that the values for the parameters in this fit are in syntony 
with the theoretical expectations for the proton impact factor since Qq is 
very similar to the confinement scale and the value of 5 sets the maximal 
contribution from the impact factor also in that region. This is reasonable 
given that the proton has a large transverse size. 



2.2 F L 

The longitudinal structure function is an interesting observable which is very 
sensitive to the gluon content of the proton. We will now present our predic- 
tions for Fl using the best values for the parameters previously obtained in the 
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Figure 2: Study of the dependence of F 2 (x,Q 2 ) on Q 2 using the LO photon 
impact factor (solid lines) and the kinematically improved one (dashed lines). 
Q 2 runs from 1.2 to 200 GeV 2 . 
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Figure 3: Study of the dependence of F 2 (x,Q 2 ) on x using the LO photon 
impact factor (solid lines) and the kinematically improved one (dashed lines). 
Q 2 runs from 1.2 to 120 GeV 2 . 
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fit of i*2. We will see that the agreement with the data is very good. First, Q 2 
is fixed and the x dependence is investigated in Fig. |4| The experimental data 



have been taken from 18 . To present the Q 2 dependence it is convenient to 
calculate, for each bin in Q 2 , the average value of x, see Fig. [5] In some sense 
this is a similar plot to the one previously presented for A in the F2 analysis 
and we can see that the effect of using different types of impact factors is to 
generate a global shift in the normalization. Again we note that we have an 
accurate description of the transition from high to low Q 2 , which was one of 
the main targets of our work. 



3 Predictions for future colliders 

While our predictions for the structure functions are in agreement with the 
data from the HERA collider experiments HI and ZEUS, these observables 
are too inclusive to provide unambiguous evidence for BFKL evolution (for 
other recent studies in this context see [Tt] ) . Comparable in quality fits can 



be obtained by both DGLAP evolution and saturation models, see e.g. 18,19 
In order to distinguish among different parton evolution pictures new collider 
experiments are needed, such as the proposed Electron-Ion- Collider (EIC) at 
BNL/JLab (USA) [20] and the Large Hadron Electron Collider (LHeC) at 
CERN (Switzerland) |2i], which will be able to measure both F2 and Fl at 
unprecedented small values of Bjorken x. In Fig. [6]we present two studies with 
our predictions for F 2 and F L down to values of x = 10~ 6 . 



4 Conclusions 

We have presented an application of the BFKL resummation program to the 
description of the x and Q 2 dependence of structure functions as extracted 
from Deep Inelastic Scattering data at HERA. We have also provided some 
predictions for these observables at future colliders. In order to obtain the cor- 
rect dependence on the virtuality of the photon at high values of the scattering 
energy, we have included in the BFKL kernel the main collinear contributions 
to all orders. We have also used optimal renormalization and an analytic run- 
ning coupling in the infrared in order to accurately describe the regions of low 
Q 2 . Our next task will be to use these parameterizations to describe more 
exclusive observables, such as heavy quark and multi-jet production, at the 
Large Hadron Collider at CERN. 
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Figure 4: Fit to Fl with the LO photon impact factor (solid lines) and the 
improved one (dashed lines). The experimental data are taken from [i~8]. 
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Figure 5: The proton structure function Fi as a function of Q 2 . The average x 
values for each Q 2 of the HI data (black) are given in Figure 13 of [18]. ZEUS 
data are taken from |19|. The solid line represents our calculation with the LO 
photon impact factor and the dashed line using the kinematically improved 
one. 



F 2 predictions for LHeC ^ predictions for LHeC 




Figure 6: Predictions for F 2 (left) and Fl (right) for LHeC. On the left plot, the 
lOGeV 2 can be compared with Figure 4.13 of 21 . Simulated 



curve with Q 2 

measurements for Fl in the kinematic range plotted here (right) can be found 
in Figure 3.7 of the same reference. 
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